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ABSTRACT 

We investigate in detail the role of active galactic nuclei (AGN) on the physical state of the gas 
in galaxy groups and clusters, and the implications for anisotropy in the cosmic microwave background 
(CMB) from Sunyacv-Zcldovich (SZ) effect. We have recently showed that AGNs can significantly change 
the entropy of the intracluster medium (ICM) and explain the observations of excess entropy in groups 
and clusters. AGNs are assumed to deposit energy via buoyant bubbles which expand as they rise in the 
cluster atmosphere and do PdV work on the ICM. Here, we include the effect of thermal conduction, 
and find that the resulting profiles of temperature and entropy are consistent with observations. Unlike 
previously proposed models, our model predicts that isentropic cores are not an inevitable consequence 
of preheating. The model also reproduces the observational trend for the density profiles to flatten in 
lower mass systems. 

We deduce the energy i?agn required to explain the entropy observations as a function of mass of groups 
and clusters Mdustcr and show that E'agn oc M^^^^^^^. with a ~ 1.5. We demonstrate that the entropy 
measurements, in conjunction with our model, can be translated into constraints on the cluster — black 
hole mass relation. The inferred relation is nonlinear and has the form Mbh oc M°[^^^^^. This scaling is 
an analog and extension of a similar relation between the black hole mass and the galactic halo mass 
that holds on smaller scales. 

In addition, we study the implications of these results for thermal SZ effect. We show that the central 
decrement of the CMB temperature is reduced due to the enhanced entropy of the ICM, and that the 
decrement predicted from the plausible range of energy input from the AGN is consistent with available 
data of SZ decrement. We also estimate the Poisson contribution to the angular power spectrum of 
the CMB from the SZ effect due to AGN heating. We show that AGN heating, combined with the 
observational constraints on entropy, leads to suppression of higher multipole moments in the power 
spectrum and we find that this effect is stronger than previously thought. The supression in the power 
spectrum in our model is due to depletion of gas from the central regions that is more efficient in low 
mass clusters and groups than in massive clusters. 



Subject headings: cosmology: theory — galaxies: clusters: general 
X-rays: galaxies: clusters 



cosmic microwave background 



1. INTRODUCTION 

The formation of structures in the Universe is believed 
to be hierarchical, as primordial density fluctuations, am- 
plified by gravity, collapse and merge to form progres- 
sively larger systems. This hierarchical development leads 
to the prediction of self-similar scaling relations between 
systems of different masses and at different epochs (Pee- 
bles, 1980). These structures contain two components - 
the gravitationally dominant dark matter and the baryons 
contained in these potential wells whose response to pro- 
cesses other than gravitational interactions bring about 
deviations from the self-similar scaling relations. 

Clusters and groups of galaxies contain dark matter and 
hot, diffuse gas called the intracluster medium. It was be- 
lieved that this intracluster gas follows a self-similar scal- 
ing relations. However, recent observations of clusters and 
groups of galaxies have shown that the scaling relations are 



not self-similar. The observed relations of different physi- 
cal parameters of the ICM such as density, temperature, X- 
ray luminosity and entropy have mostly confirmed the re- 
quirement for non-gravitational processes like AGN heat- 
ing and radiative cooling (Lloyd-Davies et al 2000, Pon- 
man et al. 2003, Sanderson et al. 2003, Pratt & Arnaud 
2003, Pratt & Arnaud 2005). Simulations and theoreti- 
cal models of clusters with gravitational processes alone 
also point to the fact that the entropy or X-ray luminosity 
observations can be matched only with non-gravitational 
heating (Roychowdhury & Nath 2003). Many theoretical 
models have been proposed to explain these X-ray obser- 
vations by heating from supernovae (Valageas & Silk 1999; 
Wu, Fabian & Nulsen 2000), radiative cooling (Bryan 
2000; Voit & Bryan 2001; Muanwong et al. 2002; Wu & 
Xue 2002a; Dave, Katz & Weinberg 2002, Tornatore et al. 
2003), accretion shocks (Tozzi & Norman 2001; Babul et 
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al. 2002), quasar outflows (Nath & Roychowdhury 2002), 
and "effervescent heating" (Roychowdhury et al. 2004; 
hereafter RRNB04). More information on various heat- 
ing models can be found in a review by Gardini & Ricker 
(2004). 

In a recent study Croston et al. (2004) presented the 
luminosity-temperature relation for groups and separated 
the sample into "radio loud" or "quiet" objects. They 
showed that "radio loud" groups deviated more from the 
self-similar scaling relation than the; radio quiet ones. ^ 
This argument adds more credibility to the idea that 
AGNs are responsible for the entropy "floor" and devi- 
ations from self-similar scaling relations. 

In a recent paper, we have studied the effect of buoyant 
bubbles which deposit energy into the ICM as they rise and 
expand in the cluster atmosphere (RRNB04). We studied 
the evolution of the density, temperature and entropy pro- 
files of the ICM under the effect of "effervescent heating" 
from buoyant bubbles, radiative cooling and convection, 
and deduced the energy required from the AGN to satisfy 
the entropy observations of the ICM. Although the con- 
straints from entropy observations at two different radii 
(0.1r2oo and rsoo) were satisfied, the temperature and en- 
tropy profiles were somewhat different from observed pro- 
flies, in that the entropy profiles had a flattened core and 
the temperature in the inner regions were enhanced (see 
RRNB04 for details). To resolve this issue, we study the 
effect of thermal conduction in this papc;r. 

The role of thermal conduction in the intracluster 
medium has generated a lot of interest in recent times. 
Its role, however, has mostly been discussed or studied in 
detail in reference to cooling flows in the central regions of 
clusters. Many authors have investigated whether thermal 
conduction alone can act as a heating source in the centre 
of clusters to stop the gas from radiatively cooling to very 
low temperatures (e.g., Zakamska & Narayan, 2003; Voigt 
et al. 2002; Loeb 2002). It has also been studied along 
with other heating mechanisms like "effervescent heating" 
again in the context of cooling flows in the centers of clus- 
ters (Ruszkowski & Begelman, 2002). In this paper, we 
study the effect of thermal conduction along with heating 
via AGN in changing the entropy profile of the ICM. 

Until recently it was only X-ray observations that 
yielded information about the entropy excess. Due to ad- 
vances in detectors and new observing strategies (Birkin- 
shaw 1999, Grego et al. 2001; Grainge et al. 2002; Reese 
et al. 2002; Zhang & Wu 2000), the thermal Sunyaev- 
Zeldovich (SZ) effect (Sunyaev & Zeldovich 1972, 1980) 
is emerging as an independent test of the density and the 
thermal structure of clusters, thus equivalently of the en- 
tropy excess. 

Many authors have investigated the role of excess en- 
tropy in clusters on the SZ effect and tried to quantify it 
(White et al. 2002; Springel et al. 2001; da Silva et al. 
2001, 2004; Cavaliere & Menci 2001: Holder & Carlstrom 
2001 & McCarthy et al. 2003a, 2003b). Holder & Carl- 
strom (2001), Cavaliere & Menci (2001) and McCarthy 
et al. (2003a) have also examined a few SZ scaling rela- 
tions for individual clusters. They have shown that the SZ 

^ The extrapolation of the cluster L-T relation onto Figure 2 of Cros 
in the radio quiet AGN could be due to the fact that they are not a 
out this possibility, see also Donahue (2005)) 



decrement is reduced in individual clusters as a result of 
energy injection and that the SZ scaling relations deviate 
from the self-similar predictions. In a more recent effort, 
Lapi et al. 2003 have estimated the enhancements in the 
SZ effect due to transient blastwaves from quasars and the 
depressions when the hydrostatic equilibrium is recovered. 

In this paper we further explore the consequences of 
heating the intracluster medium via the "effervescent heat- 
ing" mechanism (Begelman 2001, Ruszkowski & Begelman 
2002, RRNB04) and thermal conduction. We also focus on 
the SZ decrement resulting from this heating model and 
calculate the angular power spectrum of the CMB due to 
effervescent heating and thermal conduction. Using the 
entropy data, we also put constraints on the relation be- 
tween the black hole mass and the cluster mass. We show 
that it is has to be nonlinear in analogy to the black hole — 
galaxy halo relation that is expected to hold on smaller 
scales. 

The paper is organized as follows. In Section 2 we briefly 
describe our model. In Section 3 we estimate the cen- 
tral SZ decrements. We also simulate the evolution of 
the central SZ signal due to AGN heating, cooling and 
conduction. In Section 4 we estimate the angular power 
spectrum of the SZ temperature decrement in our mod- 
els. We present our results and discussion in Section 5. In 
Section 6 we discuss the relation between the mass of the 
central black hole and the mass of the cluster. Finally, our 
conclusions are summarized in Section 7. 

We assume throughout the paper = 0.71, = 
0.29, Ob = 0.047 and h = 0.71 which are the best fit 
parameters from WMAP (Spergel et al. 2003). 

2. MODEL OF THE INTRACLUSTER MEDIUM 

2.1. The default state of the ICM 

The details of the initial conditions of our model are 
similar to those adopted by RRNB04. In brief, we assume 
that the ICM is characterized by a "universal temperature 
profile" in the range 0.04 < r/rvir < 1 given by 



(T) {1+r/ay 

where (T) is the emission-weighted temperature, b = 1.33, 
a = r^ii and 6 = 1.6. 

This choice of the temperature profile as the initial tem- 
perature profile stems from the results of recent simulation 
of adiabatic evolution of clusters (Loken et al. 2002). Since 
our goal is determine the effect of energy input from AGN, 
the most natural choice for an initial profile is that which 
arises from cluster evolution without energy input. Also, 
as Roychowdhury & Nath (2003) have shown, this initial 
temperature profile results in entropy profiles closer to ob- 
servations than other profiles used previously in literature. 
The choice of this temperature profile therefore decreases 
the discrepancy in entropy and therefore provides more 
conservative estimate of the energy required from AGN. 

Density profiles are computed assuming hydrostatic 
equilibrium of the gas in the background dark matter po- 
tential. The background dark matter density profile is 
given by the self-similar Navarro, Frenk & White (NEW) 

et al. (2004) would fall above all the points. Slightly smaller deviation 
'e now but were active in the past. (We thank the referee for pointing 
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profile (Komatsu & Seljak 2001) witli a softened core (Za- 
kamska & Narayan, 2003) 

Pdra = 7 ; 77 ; To') (2) 

(r + rc)(r + Ts)^ 

where rg is the standard characteristic radius of the NFW 
profile, Tc is a core radius inside which the density profile 

is a constant and ps is the standard characteristic density 
of the usual NFW profile. The mass profile is given by 

Mdm(< r) = inpsr^w^ix), (3) 
where m{x) is a non-dimensional mass profile 



and the buoyant (effervescent) phase. It can be determined 
self-consistenly by taking into account energy losses due to 
bubble creation, 
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(L) = -r{L) 



(8) 



This comes from the fact that the energy lost within the 
radius tq due to bubble creation is PqVo where Vq is the 
volume of the bubble at the creation region or transition 
region, ro, and Pq is the pressure at that region. How- 
ever, the energy lost due to bubble expansion, i.e. in the 
effervescent phase is 



m{x) = 



(1-Xe)2 
(1 - 2Xc) 



ln(l + x/xc) 



where x = r/vg and Xc = r/vc- If Tc = 0, the usual 
mass distribution is recovered as in Komatsu & Seljak 
(2001). We follow Zakamska & Narayan (2003) and as- 
sume fc = rs/20. This is a reasonable choice as cluster 
lensing studies suggest that the core radius can be ~ tens 
of kilo-parsecs (Tyson et al., 1998; Shapiro & Ilicv, 2000). 
We investigate the effect of the smoothing of the dark mat- 
ter profile on our results. 

We do not model the effects of magnetic fields explic- 
itly. However, magnetic fields are likely to be below the 
equipartition value and thereby dynamically unimportant 
(Fabian et al. 2002). Polarization measurements suggest 
that the magnetic fields in the vicinity of bubbles have 
subequipartition values (Blanton et al. 2003). Neverthe- 
less, they may still be important for suppressing instabili- 
ties on bubblc-ICM interfaces. Such effects may implicitely 
be included in the bubble heating model. 

2.2. Heating, cooling and thermal conduction 

2.2.1. Effervescent heating and radiative cooling 

The effervescent heating mechanism is a gentle heating 
mechanism in which the cluster gas is heated by buoyant 
bubbles of relativistic plasma produced by central AGN 
(Begelman 2000, Ruszkowski & Begelman 2002). The av- 
erage volume heating rate is a function of the ICM pressure 
gradient and is given by 



(5) 



where Pgas is the ICM pressure, 7b is the adiabatic index 
of buoyant gas in the bubbles and h{r) is the normalization 
function 



Mr) = ^[l-e-'-/-k-^ 



(6) 



In equation (6), (L) is the time averaged energy injection 
rate and ro is the inner heating cut-off radius. The nor- 
malization factor q is defined by 



q = 



= Z'^"" p(7b-i)/7b 111^(1 _ e-^/^o)dr, (7) 



r d\nr 



where rmax = ^200- The inner heating cut-off radius is 
the transition region between the bubble formation region 



bubble 



dV. 



(9) 



Po 



Assuming adiabatic evolution of the gas inside the bubbles 
(low density, low radiative losses) and mass conservation in 

the bubble one has dV = (l/7b)(Po/Pbubbie)^/^'>^^o(iP/P- 
On integrating, one gets effervescent energy as 3PoVo for 
7b = 4/3. It can be easily seen from above that the en- 
ergy loss due to bubble creation is approximately 25% of 
the total energy available for heating. Thus this condition 
sets the inner cut-off radius, ro, since this 25% of the to- 
tal energy or, equivalently. 1/4(L) will be lost for bubble 
creation and injection within this radius. We find that ro 
roughly turns out to be around ~ 20 - 45 kpc depending 
on the cluster mass Mci(= Mvir) and temporal profiles of 
pressure and density. Higher mass clusters have greater 
ro- 

The volume cooling rate is calculated using a fit (Nath 
2003) to the normalized cooling function A^{T) for a 
metallicity of Z/Zq = 0.3, as calculated by Sutherland 
& Dopita (1993). This cooling function incorporates the 
effects of free-free emission and line cooling. Thus, the vol- 
ume cooling rate is F = n^A^{T), where = 0.875{p/mp) 
is the electron density. 

2.2.2. Thermal Conduction 

A potential difficulty in raising the entropy of the in- 
traclustcr medium at large radii (r > 0.1r2oo) by means 
of a central heating source is that the energy required is 
fairly large (> 10^^ erg) over a period of the age of the 
cluster. This sets up a negative entropy gradient in the 
central regions of clusters (see Figure 2 in RRNB04) and 
a rising temperature profile (see right panel of Figure 3 
in RRNB04) which are somewhat different from the ob- 
served temperature and entropy profiles. Our previous 
analysis, however, did not include the effect of thermal 
conduction which would have decreased the temperature 
gradient in the inner region and made it consistent with 
observations. We, therefore, include the effect of thermal 
conduction here to address these issues. 

Thermal conduction has been suggested to be an im- 
portant process in galaxy clusters for quite some time 
(Bertschinger & Meiksin 1986; Malyshkin 2001; Brighenti 
& Mathews 2002; Voigt et al. 2002; Fabian, Voigt & Mor- 
ris 2002). However it is not clear as to how dominant 
it would be since magctic fields would suppress the con- 
duction co-efficient by a large amount from the classi- 
cal Spitzer value. However, recent theoretical works by 
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Narayan & Mcdvcdcv (2001), Chandran & Maron (2004), 
Loeb (2002) and several others suggests that conduction 
could be as high as 10% to 20% of the Spitzer value in the 
presence of a tangled and turbulent magnetic field. Mo- 
tivated by these results we adopt the suppression factor 
/ = 0.1. 

The flux due to thermal conduction .Fcond is given by 



Fcond = -/kVT, 

where k is the Spitzer conductivity 

1.84 X lO-^T^/^ 



K 



V A 



(10) 



(11) 



with the Coulomb logarithm In A = 37, approprite for ICM 
temperature and density. 

2.3. Evolution of the ICM 

The intraclustcr gas is assumed to be in quasi- 
hydrostatic equilibrium at all times since cooling is not 
precipitous at these radii and the heating is mild. The gas 
entropy per particle is 



S = const + 



1 



7-1 



(12) 



where a = P^^Jp^^^ is the "entropy index" and 7 is the 
adiabatic index. The particle number density of the gas, 
n, is given by n = Pgg,^/ fim^. 

During each timcstep At, the entropy of a given mass 
shell changes by an amount 

AS=^K— = ^{n-r-V-F^^^,)At. (13) 
7 — 1 (J nl 

After evaluating the change in the entropy index of each 
mass shell of gas due to heating, cooling and conduction 
after a time At using equation (13), the new entropy index 
of each shell is calculated using 



a„e„(M)=ao(M) + Aa(M) 



(14) 



where c7o(M) = Pgaso/Pg^o is the default entropy index. 
The system relaxes to a new state of hydrostatic equilib- 
rium with a new density and temperature profile. After 
updating the function a{M) for each mass shell, we solve 
the equations 



dP,^^ _ GM{< r) 
dM 

dr 

dM 



1 fP,UM)V'^''^ 



a{M) 



(15) 

(16) 



to determine the new density and temperature profiles at 
time t + At. The boundary conditions imposed on these 
equations are that (1) the pressure at the boundary of the 
cluster, r„„t, is constant and is equal to the initial pres- 
sure at Taooj i.e., -P(?'out) — constant — -Pgaso(^2oo), and 
(2) the gas mass within rout at all times is the mass con- 
tained within Tsoo for the default profile at the initial time, 
i.e., M,(r,„,) = M,„(r,oo) = 0.1333Md„(r2„„). It is impor- 
tant to note here that r^^t increases as the cluster gas gets 
heated and spreads out. 

The observed gas entropy <S(r) at O.lraoo and at r^oo is 
then calculated using 



S{r)^T{r)/nl'\r). 



(17) 



The updated values of a{M) and pressure of the ICM 
P%B.„{i') are used to calculate the heating and cooling rates 
and the conduction flux for the next time step. This is 
continued for a duration of fhoat- After that the heating 
source is switched off, putting H = 0. The cooling rate 
and conduction flux continue to be calculated to update 
the function cr(M) at subsequent timesteps, and the hy- 
drostatic structure is correspondingly evolved for a dura- 
tion of — thcat. where t^ is the Hubble time. Note that 
Tout decreases during this time since the intracluster gas 
loses entropy and shrinks. The only free parameters in 
our calculation are the energy injection rate (i) and the 
time ihcat over which the "effervescent heating" of the ICM 
takes place. After evolving the gas for the total available 
time, ~ 1-35 x 10^" years, we check whether the en- 
tropy at O.lraoo and rr^a match the observed values, and 
adjust parameters accordingly. In this way, we explore the 
parameter space of {L) and ti,oat or rather a single free pa- 
rameter, i.e. the total energy, i?agn = {L)theat for different 
cluster masses so that the entropy (after 1.35x10^° years) 
at O.lraoo and Vaoo matches the observed values. 

For numerical stability of the code, the conduction term 
is integrated using timesteps that satisfy the appropriate 
Courant condition. The Courant condition for conduction 



IS 



At, 



< 0.5 



(Ar)^nfcb 
k(7 - 1) ■ 



(18) 



The timesteps. At, used in equation (12) to update the 
entropy of the gas and calculate its pressure, temperature 
and density profiles always obey the above Courant con- 
dition (Ruszkowski & Begelman 2002; Stone, Pringle & 
Begelman 1999). 

3. THERMAL SUNYAEV-ZELDOVICH EFFECT 

The temperature decrement of CMB due to the SZ ef- 
fect is directly proportional to the Compton parameter 
(y). For a spherically symmetric cluster, the Compton 
parameter is given by 



y 



2^ / p,{r)dl 



(19) 



where ctt is the Thomson cross-section, and pdr) = 
ne(r)fcbTe(r) is the electron pressure of the ICM, where 
ne(r) = 0.875(pga,s/TOp) is the electron number density, 
fcb is the Boltzmann constant, and Tc(r) is the electron 
temperature. The integral is performed along the line-of 
-sight (l) through the cluster and the upper limit of the in- 
tegral (+i?) is the extent of the cluster along any particular 
line-of- sight. We do not include the effects of beam size in 
calculating the y parameter. This approximation is justi- 
fied by the fact that the pressure profiles are relatively flat 
in the inner region. The variation of pressure integrated 
along the line-of-sight as a function of the projected ra- 
dius is even flatter thus providing more justification for 
the above approximation. 

The angular temperature profile projected on the sky 
due to SZ effect, AT{6)/Tcmb is given in terms of the 
Compton parameter in equation (19) 
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— ^ = g{x)y{e), (20) 

JCMB 

where g{x) = a;coth(x/2)-4, x = hv/k^TcMB., 3cmb = 
2.728 (Fixsen et al. 1996). In the Rayleigh- Jeans ap- 
proximation, g{x) « —2. Wc only cvahiatc "central" SZ 
decrement from the pressure profiles of our models. In this 
case, the integral in equation (19) reduces to 

2/0 = 2^ / p,{r)dr (21) 

In the Rayleigh-Jeans part of the CMB spectrum, the devi- 
ation from the black-body spectrum results in a decrement 
of the CMB temperature, 

AT^wo«-5.5yif (22) 

We use the pressure profiles resulting from our model to 
calculate the central SZ decrement in the temperature of 
the CMB. 

4. ANGULAR POWER SPECTRUM 

The angular two-point correlation function of the SZ 
temperature distribution in the sky is conventionally ex- 
panded into the Legendre polynomials: 




(23) 

Since we consider discrete sources, we can write Ci = 

(p) (c) (P) 
C\ +C\ , where ' is the contribution from the Pois- 

son noise and is the correlation among clusters (Pee- 
bles 1980, § 41). We define the frequency independent part 

in the power spectrum as C**-''^'' = Cijg^ix). The integral 

expression of ' can be derived following Cole & Kaiser 
(1988) as 

(24) 

where V{z) is the co-moving volume and yi is the angular 
Fourier transform of y{6) given by 

yi = 2Tr J y{0)M{l + l/2)e]ede, (25) 

where Jo is the Bcssel function of the first kind of the inte- 
gral order 0. In equation (24), z^cc is the redshift of pho- 
ton decoupling and dn/dM is the mass function of clusters 
which is computed in the Press-Schechter formalism (Press 
& Schechter 1974). The mass function has been computed 
using the power spectrum for a ACDM model with nor- 
malization of as = 0.9. We choose Mmin = 5 x lO^^M© 
and Mmax = 2 x IO^^Mq and integrate till redshift of 
z = 5 instead of Zdcc- This is done because the integral in 
equation (24) is found to be insensitive to the upper limit 
in redshift beyond z = 4, the reason being that the mass 
function is exponentially suppressed beyond that redshift 
in this mass range. 



5. RESULTS 

In this section, we discuss our results for cluster evo- 
lution due to heating, cooling and conduction. We also 
discuss our results for the central SZ decrement for clus- 
ters with masses ranging from Md = 5 x 10-^^-2 x lO^^M©. 

The gas is heated for a time thcat and cooled simultane- 
ously. After this time, the heating source is switched off. 
The gas is then allowed to cool radiatively until a total 
simulation time of tH=l-35xlO^° years has elapsed. The 
final entropy values at O.lrjoo and rgoo are compared with 
the observed ones. 

In Figure (1), the evolution of the density and temper- 
ature profiles of the ICM are shown for a cluster of mass 
6 X 10" Me and for a luminosity of (L) = 5.25 x 10^^ 
ergs~^. The gas density decreases with time during the 
heating epoch, and increases due to radiative cooling and 
conduction after the heating source is switched off. It is 
interesting to note that the changes in density are mini- 
mal beyond 0.5r2oo, as compared to 0.2r2oo in Figure (3) in 
RRBN04, and that conduction plays a very important role 
in regulating the density profiles after the heating source is 
switched off. It is seen that conduction actually decreases 
the density of the gas at larger radii (beyond 0.5r2oo) by 
conducting heat out from the central regions. This is seen 
more clearly if one studies the evolution of the tempera- 
ture profiles. After the heating source is switched off, it is 
seen that the temperature of the central regions fall very 
rapidly since conduction pumps out heat from the cen- 
tral regions and redistributes it in the outer regions of the 
cluster. Thus the temperature profiles do not rise towards 
the centre as compared to what is seen in Figure (3) in 
RRBN04. On the other hand, their evolution shows a rise 
in the outer regions (beyond 0.5r2oo) due to thermal con- 
duction even after the heat source has been switched off. 
Thus conduction acts like a heating source for larger radii. 

Figure (2) shows the time evolution of scaled entropy 
profiles of a cluster of mass M^, = 6 x 10^'* Mq for (L) = 
5.25x10^^ erg s""'^. We use the same method of emissivity 
weighting as in Roychowdhury & Natli (2003) to calculate 
the average quantities. The entropy profiles are plotted in 
time-steps of 5 x 10^ years. They are seen to rise with time 
as the ICM is heated. Then, after the heating is switched 
off (after thoat = 5 x 10^ years), the gas loses entropy due 
to cooling and the profiles are seen to fall progressively. 
The inclusion of conduction removes the negative gradient 
in the scaled entropy profiles in the central regions of the 
cluster (within 0.5r2oo- These entropy profiles do not show 
any flat entropy core unlike in RRBN04 (left panel in their 
Figure 2). Thus, this probably indicates that thermal con- 
duction is a more plausible process in ICM than covection 
for such gentle AGN heating. 

The model is constrained by entropy-temperature rela- 
tion. At fixed cluster temperature this corresponds to a 
given density (at both radii for which the entropy data is 
provided). Now, since at a given density we have an ad- 
ditional constraint from the L-T relation (that our model 
fits reasonably well; see below), we iniplicitely satisfy the 
constraints on the slope of the density profile. Thus, as the 
model fits both the entropy data and the observed L-T re- 
lation that specify the slopes of the density profiles, these 
slopes must also be consistent with observations that show 
flattening in low mass systems. Indeed, this flattening is 
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Scaled Radius {r/r^J Scaled Radius {r/r^J 

Fig. 1. Gas density (left panel) and temperature (right panel) profiles as a function of scaled radius (r/r2oo), for the cluster 
model described in the text. Dashed lines represent density profiles when heatiug is active and the dot-dashed line represents 
density profiles after the heating source has been switched off (i.e. after theat > 5 x lO'' years). The profiles are plotted after 
every 5 x 10* years till Hubble time. It is seen here that conduction removes temperature gradients in the central regions (within 
0.2r2oo) and flattens the temperature profile. This is the reason thermal conduction was included in the cluster model as compared 
to RRNB04, wherein left panel of Figure (3) shows rising temperatures in the centre. The density profiles are seen to rise after the 
heating source is switched off and thermal conduction and radiative cooling are the only two processes which are active. Initial 
density and temperature profiles correspond to the solid curves in both plots. 




0.1 1 

Scaled Radius (r/rg^g) 

Fig. 2. — Scaled entropy profiles as a function of scaled radius for a cluster of mass 6 X 10^* Mq heated by an AGN with (L) = 3 X 10''^ 
erg s~^. The scaled entropy profiles are plotted at intervals of 5x10® years. The dashed lines correspond to times when heating is on and 
the dot-dashed lines correspond to times when heating has been switched off. They are seen to rise as the gas is heated and then fall as the 
gas cools. It is seen here that the final entropy profiles (after Imbble time) has neither negative entropy gradient nor any entropy core in the 
central regions, as compared to Figure (2) in RRNB04. Initial states correspond to the solid curves. 
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Fig. 3.— Final density profiles for objects of mass 4.5 x lO^^ M0 (solid line), 2.0 x lO^* Mq (short dashed line), 9.0 x lO" Mq (long 
dashed line). The assumed heating time theat ™ this example was 5 x 10^ years. Note the flattening in the density profile as the mass of the 
system is decreased. 
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Fig. 4. — This figure shows the permitted total injected energy range as a function of the cluster mass for ICM heating times between 
iheat = 5 X lO^yr (upper envelope) and t/ieat = 5 X lO^yr (lower envelope). The shaded region corresponds to values of i,agn that are able to 
match the entropy observations at both 0.1r2oo and rsoo. The thin solid line represents a non-linear relation between the total energy injected 
into the cluster by AGN and the mass of the cluster with an exponent of 1.5. The thick long-dashed line represents a non-linear relation 
between the total energy injected into the cluster by AGN and the mass of the cluster with an exponent of 5/3. The permitted parameter 
space comes from the sum of permitted regions that satisfy the entropy constraints at both radii for fixed t/jeat- 
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apparent in Figure 3 that shows final density profiles for 
different masses. 

We now discuss the permitted range in the total energy 
injected into the cluster, E^g„ = {L) x theat, required to 
match the observed entropy as a function of the cluster 
mass. 

Figure (4) shows the permitted total injected energy 
range as a function of the mass of cluster for heating times 
between theat = 5 x 10^ years and theat = 5 x 10^ years. 
Here the entropy is required to match obscTvations at both 
O.lraoo and r^oo- The thick solid line represents a relation 
between the total energy injected to the cluster by AGN 
and the mass of the cluster (see next section for more de- 
tails). 

Figure (5) shows the X-ray luminosity {Lx versus 
emission-weighted temperature (Tx)relation in clusters. 
The data points have been compiled from Arnaud & 
Evrard (1999) (represented by stars), Markcvitch ct 
al. 1998 (represented by open squares) and Helsdon & Pon- 
man 2000 (with error bars). The X-ray luminosity has 
been calculated within the cluster volume of 0.3r2O05 as 
done for the data. It is also seen that the X-ray luminosity 
docs not change much (within 1%) if the volume increased 
from 0.3r2oo to r2oo- The shaded region in the plot corre- 
sponds to the predicted X-ray luminosity when the cluster 
is heated by an AGN for 5 x 10^ < theat < 5 x 10^ years 
with luminosities which correspond to the two bonding 
lines of the shaded region in Figure (3). It is seen that 
the predicted luminosities of the heating model satisfy the 
data points in the low mass end as well as the high mass 
end. The solid line in the plot shows the predicted lumi- 
nosity due to the universal temperature profile and the de- 
fault density profile. We note that the X-ray luminosity is 
over-predicted by the universal temperature profile which 
indicates that the addition of non-gravitational heating is 
required to lower the X-ray luminosity to satisfy the data 
points (see Roychowdhury & Nath 2003, for more details). 

In Figure (6), the central SZ temperature decrement 
AT^wO is plotted as a function of the emission-weighted 
temperature of the cluster (T). The data points are a 
compilation of data sets from Zhang & Wu (2000) and 
McCarthy et al. (2003b). The solid line shows the pre- 
dicted ATjiiwO from the default temperature profile and 
NFW potential. The dot-dashed line shows the predicted 
AT^wo for the same temperature profile but for smoothed 
NFW potential with Tc = rs/20. The dashed line shows 
the prediction from the self-similar profile (Wu & Xue, 
2002b; Bryan 2000). 

It has been shown by Roychowdhury & Nath (2003) 
that the density profile of gas is much flatter in compar- 
ison to the self-similar profile when it assumes the "uni- 
versal temperature profile" and the standard NFW profile 
is assumed. As a result, the predicted central tempera- 
ture decrement (AT^wo) is lower than that predicted by 
the self-similar model. The normalization of AT^wO for a 
smoothed NFW profile with a core radius Tc = rs/20 is 
even lower. This happens because the introduction of a 
core radius in the dark matter profile makes the ICM den- 
sity profile shallower in the central regions as compared 
to the ICM density with a standard NFW profile. These 
decrements are closer to the data than predicted by earlier 
self-similar models for rich clusters. 



Next we examine the effects of the effervescent heating, 
radiative cooling and conduction on the central SZ decre- 
ment. We evaluate AT^^o for clusters in our sample after 
they have been evolved for a Hubble time in- The heat- 
ing source was active for theat ^h- The values of (L) 
and theat have been chosen so as to satisfy observational 
constraints on ICM entropy after tjj at 0.1r2oo and r^oo 
(Ponman et al. 2003; see Figure 4 and 5 in RRNB04). In 
other words, there is a range of (L) that satisfies the en- 
tropy observations at a l-o" uncertainty level that we used 
in our calculations. We used smoothed NFW profile with 
the core radius rc = rs/20. 

In Figure 7, the shaded region is the expectation for 
the SZ central decrement when the gas is heated by the 
central AGN. The shaded region represents the spread in 
(L) which satisfies the entropy requirements at both radii, 
0.1r2oo and rsoo for 5 x 10* < theat < 5 x 10^ years. The 
solid line shows the prediction from self-similar profile (Wu 
& Xue, 2002b; Bryan 2000). 

In Figure 8, the evolution of the central SZ decrement 
^I)iwO is shown as a result of heating, cooling and con- 
duction for a cluster of mass Md = 6 x 10^^ Mq. The 
dashed line is the result of heating for theat = 5 x 10* 
years and the solid line is the result of heating the ICM 
for theat = 5 X 10^ years. It can be seen that, as long as 
the heating source is active, AT^^o decreases. When the 
source is switched off, AT^wo starts to increase as the gas 
evolves only due to radiative cooling and conduction. This 
happens because the density, or equivalently the electron 
pressure, decreases when the gas is heated but becomes 
larger when it is allowed to cool. In addition, we have 
also plotted the default AT^^o for the same cluster with a 
point denoted with an open circle. The central SZ decre- 
ment corresponding to our heating model is always lower 
than the default value. 

Finally, we evaluate the Poisson contribution to the an- 
gular power spectrum of the SZ and compare our model 
predictions with earlier ones from self-similar models. In 
Figure (9), the thin solid line represents the angular power 
spectrum (Poisson) for the universal temperature profile 
and the corresponding density profile (Roychowdhury & 
Nath, 2003). The dashed line is for self-similar model (Ko- 
matsu & Kitayama, 1999). The shaded regions represent 
the angular power spectra from heating model. The region 
shaded with gray represents the result due to heating for 
theat = 5 X 10^ years and the region shaded with dashed 
lines represents the result due to heating for theat = 5 x 10* 
years. The spread reflects the fact that there is a range of 
(L) for a given theat that satisfies the observational entropy 
data at 0.1r2oo and ^500. 

6. THE HALO-BLACK HOLE MASS RELATION 

The relation between the mass of the group or cluster 
halo and the total injected encgy (Figure 4) can be trans- 
lated to the halo-black hole mass relation. We follow the 
arguments in Wyithe & Loeb (2003) to derive this relation. 

Mass of a virialized halo can be expressed in terms of 
its circular velocity Vc (Barkana & Loeb 2001) 

Mhaio = 4.7 x 10^^ f A Mq (26) 

V lO'^ km s / 

We can combine the above relation with the self-regulation 
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T (keV) 

Fig. 5. — Relation between bolometric X-ray luminosity Lx and emission- weighted temperature ((T)). The data points represented by 
'stars' show measurements of clusters with insignificant cooling flows compiled by Arnaud & Evrard (1999). Open squares show cooling 
flow-corrected measurements by Markevitch ct al. (1998). The data points with error bars show group data from Helsdon & Ponman (2000). 
The shaded region represents X-ray luminosity calculated for the region of Eagn shown in Figure (4) which satisfies the observed entropy 
requirements at both radii for 5 X 10** < theat < 5 X 10'^ years. The solid line represents the X-ray luminosity calculated using the default 
model of the ICM ie. the universal temperature profile. The models assume a ACDM cosmology with Q,m = 0.29, Ha = 0.71, and Oj, = 
0.047, and a Hubble parameter of ft. = 0.71 has been applied to the models and the data. 




Emission — weighted temperature (T^ i<eV) 

Fig. 6. — Observed and predicted A7)iw0 -(T) relation of clusters. The solid line represents the predicted AT^wO with "universal" temper- 
ature profile (Loken ct al. 2002), dark matter density profile given by NFW with rc = and the resulting density profile (Roychowdhury & 
Nath 2003). The dash-dotted line is the result of the dark matter profile given by NFW with rc = r's/20 and the ICM temperature profile as 
before. The dashed line is the result of self-similar profile (Wu & Xue 2002b). The data points are from Zhang & Wu (2000) and McCarthy 
et al. (2003b). 
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Emission — weighted temperature (T^ i<eV) 

Fig. 7. — Observed and predicted AT,j^v()^(7") relation of clusters. The solid line shows AT^^o for self-similar profile (Wu & Xue 2002b). 
The shaded regions are the SZ decrciuciit due to heating of the ICM, the region shaded with dashed line representing the spread in (L) which 
satisfy the entropy observations at both radii, 0.1r2oo and rsoo for theat = 5 x 10^ years and the gray region representing the spread in (L) 
for tfieat = 5 X 10^ years. The data points are from Zhang & Wu (2000) and McCarthy et al. (2003b). 
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Fig. 8. — The evolution of the central SZ temperature decrement, AT^„Oi is shown for two different Lagn corresponding to the two values 
of tfteat for A^ci = 6 X 10^'* M0. The point denoted by an open circle is the value of AT^wO for the default profile with smoothed NEW ie. 
with a core for the same mass cluster. 
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I 

Fig. 9. — The Poisson contribution to the angular power spectrum (C*) of the SZ fluctuations is plotted here as a function of /. The 
thick solid line is for the primary temperature anisotropy expected in the Rayleigh-Jeans band. The thin solid line is for the angular power 
spectrum (Poisson) due to the "universal temperature profile" of the cluster ICM and the dashed line is the expectation from a /3— profile 
and isothermal temperature profile. The dotted and the dot-dashed line represent the heated ICM models. The dot-dashed line being for 
iheat = 5 X 10* years and the dotted line is for tjieat = 5 X 10^ years. 



condition (Wyithe & Loeb 2003) 



0.5(Ob/n^)Mhaiot;g 

^dyn 



(27) 



where L-^dd is the Eddington luminosity of the central 
black hole, r] is its Eddington fraction, Fg is the fraction 

of energy generated by the black hole that is deposited 
in the gas, Qb and Qm are the mass density in baryons 
and in mass relative to the critical density, respectively. 
We assumed that the dynamical time of the halo gas is 



tdy 



r/vc, where r^ir is the virial radius of the col- 



lapsed halo (Barkana & Loeb 2001). Combining equations 
(25) and (26) we get 



Mbh = 2.2 X 10* 



VO.05 



0) 



(28) 



where M14 = Mhaio/IO^^M©. Note that the above rela- 
tion between the mass of the cluster and the central black 
hole is nonlinear. It is interesting that a similar scaling 
relation would be expected in the case of the galaxy halo 
mass Mgh and black hole mass relation. In such a case, 
the mass of the black hole scales like Mbh oc w", where 
4 $ a ^ 5. As oc Mgh/rgh and rgh oc Mgh, we have 
Mbh oc M"J^ . In fact, based of the M — a relation and 

©n 

cosmological simulations, Ferrarese & Ford (2005) derive a 
relation between the mass of the black hole and the mass 
of the galactic halo of the form M cx M^^^ . The slope of 
this relation is very similar to the one considered here. 

Asstiming that that a fraction e of black hole mass is 
converted to energy, we can also get the relation between 
the injected energy Eagn and the cluster mass 



^;a.n = 6xl0-(^) 



yFg 

0.05 



M. 



5/3 
14 • 



(29) 



This relation is the same as the one denoted by thick 
dashed line in Figure 4. We stress that the scaling of 
black hole mass and the mass of the cluster derived here is 
not unique. As Figure 4 demonstrates, relations that have 
slightly different slope are also acceptable. For example a 
somewhat shallower slope of 1.5 is also consistent with the 
entropy constraints. Nevertheless, we emphasize that suc- 
cessful fit requires the relation to be nonlinear with index 
steeper than unity. 

7. DISCUSSION 

In this paper, we have examined the effects of efferves- 
cent heating by AGN in clusters with thermal conduction 
and cooling in the context of the excess entropy require- 
ments at large radii. We have also examined the con- 
sequences of this heating, cooling and conduction on SZ 
temperature decrement. 

As is clear from Figure (4) , it is possible to heat the IGM 
with a single central AGN to match the entropy require- 



ments at both O.lronn and 



However, in order to match 



the entropy at both radii, the total injected energy ii^ag„, 
for a given value of theat ^ tn, must be tightly constrained. 
In fact, our calculations have shown that for any value of 
ihoat < ^H, i.e., for any heating time (or AGN lifetime), it 
is always possible to satisfy the entropy observations at 
both radii with a single value of the luminosity (L) . This 
is different from the cooling flow problem, where (L) must 
be finely tuned to match the cooling rate (Ruszkowski & 
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Begelman 2002), because cooling effects on large scales are 
rather mild and, thus, the results depend mostly on the 
total injected energy E^^„ = (L) x theat- Thus, if we can 
fit the observed entropy values for just one pair (L) and 
theat, we can do so for a wide range of such pairs. 

We note here that the inclusion of thermal conduction 
brings down the energy which has to be provided by the 
AGN over its life-time to satisfy the entropy observations 
at both radii for all heating times as compared to our 
model in RRNB04. In addition, for shorter heating times, 
the Eagn is less or comaparable to the energy pumped in 
for longer heating times. This is in contradiction to our 
findings in RRNBG4. This happens here because ther- 
mal conduction acts as a heating source after the AGN 
is switched off (for tn — theat) and raises the entropy at 
large radii (at r = rsoo)- The results are mostly sensitive 
to the total energy input from the black hole, rather than 
to (L) and theat separately. As a consequence, satisfactory 
fits can be obtained as long as the total injected energy 
falls within a relatively narrow range of values, which de- 
pends on the cluster mass (Figure (3)). 

Finally, we note that cooling and thermal conduction 
play important roles in controlling the heating mechanism 
so that the entropy profiles broadly match the observed 
entropy profiles in clusters (Ponman et al. 2003). Notably, 
in the later stages of evolution of the gas, after the heat- 
ing source is switched off, conduction removes negative en- 
tropy gradients in the central regions of the cluster. More- 
over, there is no entropy core seen in the final stages of the 
evolution of the ICM. Instead we see positive entropy gra- 
dients as observed in the entropy profiles of galaxy groups 
(Mushotzky et al. 2003, Ponman et al. 2003). Unlike pre- 
viously proposed models, our model predicts that isen- 
tropic cores are not an inevitable consequence of preheat- 
ing. However, the clusters that show isentropic core have 
also been observed (Ponman et al. 2003). We note that 
our entropy profiles show a core while the source of heat- 
ing is active. It is conceivable that the clusters which show 
entropy cores are being observed during the active phase 
of the AGN duty cycle. 

As seen in Figure (6), the central SZ decrement for indi- 
vidual clusters is diminished if the "universal temperature 
profile" corresponding to the gravitational interactions of 
the ICM with the background dark matter is used. This 
reduces the discrepancy between the predicted central SZ 
signal and the observations compared to the predictions 
based on the standard self-similar profiles. The solid line 
in Figure (6) matches the observed points for rich clus- 
ters ((T) < 5 keV) better than the dashed line, indicating 
that the requirement of non-gravitational heating for rich 
clusters is lower than previously thought. It is important 
to note here that the introduction of a core radius in the 
dark matter profile brings the SZ temperature decrements 
down further to match the observed points. This happens 
because the pressure profile of the intracluster medium be- 
comes shallower in the central regions than for the NFW 
dark matter profile without smoothing. 

We also find that the inclusion of effervescent heating, 
cooling and conduction have a significant effect on the SZ 
signal. The gas in the central region is depleted as a re- 
sult of effervescent heating and thus we see a diminution 
in the SZ signal as a result of AGN heating. Similar con- 



clusion was reached by McCarthy et al (2003a, 2003b) and 
Cavaliere & Menci (2001). However, in our case, there is 
a spread in the SZ signal. The spread reflects the fact 
that there is a range of energies required to satisfy the ob- 
served entropy at O.I-R200 and -R500 for a particular theat- 
Shaded regions in figure 2 are based on fits of the heat- 
ing model to the observational entropy data. The region 
shaded with dashed lines corresponds to the predicted SZ 
signal for the range in ii^agn that satisfies the entropy re- 
quirements at both radii for theat = 5 x 10*^ years and 
the region shaded gray corresponds to the range in i^agn 
that satisfies the entropy requirements at both radii for 
theat = 5 X 10^ years. As seen in Figure 7 the data for 
SZ are still not good enough to constrain heating models 
at present. We hope that future SZ observations of low 
temperature clusters will constrain the models better in 
the regime where discrepancies between self-similar and 
heating models are more pronoimced. The capabilities of 
future experiments like SZA are such that a decrement in 
the CMB temperature of 10 /zK could be detected (S. Ma- 
jumdar, private communication). This would mean that 
the lowest mass systems that SZA could observe would 
have T~ 2 kcV. It should be noted here that the parame- 
ter Lagn and thus i?agn estimated in this paper are slightly 
different from those in RRNB04. This is due to differ- 
ences in the effervescent heating mechanism and the dark 
matter profile. The inclusion of a self-consistent method 
for estimating the inner heating cut-off radius, ro and a 
smoothed NFW profile with a core radius have changed 
the range in E^gn that satisfies the entropy requirements 
at two different radii, 0.1r2oo and r^oo, for the two values 
of theat stated above. 

Wc have also studied the time evolution of the SZ tem- 
perature decrement due to AGN heating, cooling and con- 
duction. It is seen in Figure (8) that the SZ signal is 
diminished as a result of AGN heating at all times in com- 
parison to the default SZ signal. This is in contrast to the 
transient phases with enhancement of SZ signals predicted 
by Lapi et al. (2003) owing to strong feedback mecha- 
nisms that they assumed. Also, in the case of effervescent 
heating, the SZ decrement is lower than the default case 
even after the gas has evolved for a time tn — theat after 
the heating has been switched off. We note here that the 
values of (L) and theat have been chosen so as to satisfy 
the entropy requirements after evolving for tn- 

Finally, we examine the effect of the universal temper- 
ature profile and AGN heating on the Poisson part of the 
angular power spectrum. We note that the peak of the SZ 
power spectrum is somewhat sensitive to the amount of 
heating that is added to the cluster gas. The peak of Ce is 
at a lower £ for a larger theat and at a higher i for a smaller 
theat- Also, the effect due to heating is larger for higher 
values of i. In our case the suppression takes place because 
the gas depletion from the central regions is more efficient 
in low mass groups than in rich clusters. Therefore, the SZ 
signal is suppressed more efficiently at smaller scales and, 
thus, larger i. We note that the general trend for the Pois- 
son contribution to the spectrum to be supressed at higher 
£ has been noted by others (Holder & Carlstrom, 1999; 
Kitayama & Komatsu 1999; Springel, White & Hernquist, 
2001; Holder & Carlstrom 2001, Zhang & Wu 2003). This 
effect leads to a mismatch between the observed and the- 
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oretical spectra (Dawson ct al. 2001, Mason et al. 2003) 
when preheating required to account for the entropy floor 
is considered. In our physicaUy motivated effervescent 
AGN heating model, that fits the X-ray observations of 
entropy of the ICM, the power spactrum at small scales 
is even lower than previously thought. This suggests that 
other sources, such as, e.g., "dead" radio galaxy cocoons at 
higher redshifts (Yamada, Sugiyama & Silk 1999), should 
significantly contribute to anisotropics in the cosmic mi- 
crowave background at large We emphasize that our 
model deals only with global, average effect of heating on 
clusters and neglects small scale fluctuations in the gas 
distribution that are associated with heating. 

8. CONCLUSION 

The primary aim of this work was to study the implica- 
tions AGN heating and thermal conduction on the global 
properties of groups and clusters of galaxies including the 
Sunyaev-Zeldovich effect. 

We have demonstrated that the available entropy data, 
in conjuction with our feedback model, put constraints on 
the relation between the total energy injected by the AGN 
and the mass of the cluster (Sagn oc M^^^^^^). The inferred 
black hole-halo mass scaling {Mbh oc M^^^j^j., a ^ 1.5) is 
an analog and extension of the similar relation between 
the black hole mass and the mass of the galaxy halo that 
holds on smaller scales. 

The "effervescent heating" mechanism heats the gas in 
the central regions of clusters and makes the gas density 
profile shallower. This reduces the electron pressure of 
the gas and, thus, reduces the SZ temperature decrement. 
This is in accordance with the findings of other authors 
(Cavaliere & Menci 2001, McCarthy et al. 2003b, Lapi at 
al. 2003). Here we have also shown how the SZ decrement 
would evolve as heating, cooling and conduction regulate 
the physical state of the ICM. Our heating model is con- 
sistent with the available entropy data at 0.1r2oo and r5oo. 
We give specific predictions of our model for the SZ decre- 
ment for low mass (low temperature) clusters. Future ob- 
servations performed with, e.g., Sunyaev-Zeldovich Array 
(SZA) or Combined Array for Research in Millimiter-wave 
Astronomy (CARMA) will be able to test these predic- 
tions. 

We also point out that the "universal temperature pro- 
file", that takes into account pure gravitational interac- 
tions, leads to lower SZ decrement than that calculated 
assuming that the ICM has a self-similar profile and is in 
better agreement with data for high (T) clusters. This im- 
plies that the discrepancy between observations and mod- 



els without heating is reduced. 

Unlike in the case of previously proposed models, we 
found that isentropic cores are not an inevitable conse- 
quence of preheating. This is consistent with observations 
of groups that do not show large isentropic cores (Pon- 
man et al. 2003, Mushotzky 2003). Clusters that show 
isentropic core have also been observed (e.g., Ponman et 
al. 2003). We note that our entropy profiles show a core 
while the source of heating is active which may explain 
suc;li casc;s as well. It is conceivable that the clusters which 
show entropy cores are being observed during the active 
phase of the AGN duty cycle. This suggests that there 
may exist an observationally testable correlation between 
the presence of isentropic cores and radio emission from 
the central galaxy. However, we note that such a correla- 
tion may not be strong as the radio emission may quickly 
fade away while the mechanical effects of feedback may 
still be present. 

We have also demonstrated that the model reproduces 
the observed trend for the density profiles to flatten in low 
mass systems. 

Finally, we also estimated the angular power spectrum 
of the CMB due to the SZ effect from Poisson distributed 
clusters. We showed that the average effect of heating 
is to reduce the SZ signal and thus the angular power 
spectrum. The flnding that the power spectrum at large 
/ is suppressed is consistent with previous results (e.g., 
Komatsu & Kitayama 1999, Holder & Carlstrom 2001). 
However, our results indicate that the contribution to the 
power spectrum that results from AGN heating, that is 
consistent with entropy measurements, is lower than pre- 
viously thought. This suggests that other sources, such as, 
e.g., dead radio galaxy cocoons at higher redshifts (Ya- 
mada, Sugiyama & Silk 1999), should significantly con- 
tribute to small scale anisotropics of the cosmic microwave 
background. 
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